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The effects of lattice relaxation in LaTiOa/SrTiOs superlattices are investigated using a combina- 
tion of LDA+t/ density functional theory, and Hartree-Fock effective Hamiltonian calculations. We 
find noticeable (~ 0.1-0.2 A) distortions of the TiOe octahedra in the near-La region. The resulting 
screening changes the Ti d-electron density substantially. Tight-binding fits to the relaxed-lattice 
band structure, combined with Hartree-Fock calculations of the resulting model, reveal a novel 
phase with xy orbital order, which does not occur in bulk LaTi03, or in the hypothetical unrelaxed 
structure. 
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Recent advances in the techniques of pulsed laser de- 
position and molecular beam epitaxy have allowed the 
creation of "oxide heterostructures" consisting of alter- 
nating layers (of arbitrary number of unit cells) of dif- 
ferent transition metal oxide compounds. The physics 
becomes particularly interesting when one or more of the 
constituent layers is a compound with correlated electron 
properties 0,12 such as Mott insulating behavior, high 
temperature superconductivity or colossal magnetoresis- 
tance. In addition to advances in fabrication, recent years 
have seen tremendous p rog ress in the characterization 

Hil Il3| and theoretical analysis 
of such structures. However, be- 
fore this work, the crucial issue of lattice relaxation has 
not been systematically addressed. 

We expect the lattice relaxation issue to be important 
because a key feature of heterostructures is "charge re- 
distribution" 0; that is the redistribution of electrons 
from one layer to another, and consequent changes from 
the bulk valence values, driven by the difference in elec- 
trochemical potentials between the component materi- 
als. Simple estimates and tight-binding model studies 
suggest that the resulting dipole layer involves a charge 
density of about he / unit cell, displaced by one to two 
lattice constants [l4j. The electric field associated with 
such a dipole layer is not small, ~ 0.1 eV/A, and may be 
expected to drive significant changes in atomic positions 
compared with the bulk materials. Further, many exper- 
imentally relevant heterostructures involve ferroelectric 
or nearly ferroelectric materials, for example SrTiC>3, for 
which enhanced dielectric response is an issue. 

The purpose of this paper is to calculate the magnitude 
and nature of the lattice relaxation in LaTi03/SrTiC>3 
heterostructures, and to determine its effect on the 
electronic properties. First, we calculate the relaxed 
structures of [LaTi03]„[SrTi03] m superlattices using the 
LDA+t/ method of density functional theory within the 
projector augmented wave (PAW) approach 2Qj as im- 



plemented in the Vienna Ab initio Simulation Package 

(VASP) [Hill- We consider I 001 ] ( n " m ) superlattices 
in which a unit consisting of m planes of LaTiC>3 fol- 
lowed by n planes of SrTiC>3 layers is repeated in the [001] 
(z) direction, similar to those studied experimentally by 
Ohtomo et al. . Technical details include the rotation- 
ally invariant LDA+t/ method of Liechtenstein et al. [3] 
with U = 5 and J = 0.64 eV for the Ti d states Hi]. In 
addition, if we treat the La / states within the LDA, we 
find that the empty La / bands lie only ~ 2 eV above 
the Fermi level, leading to a spurious mixing and level 
repulsion with the Ti-derived d bands. Since in practice, 
the La / bands should lie much higher in energy [25| , we 
impose a large U of 11 eV, and J = 0.68 eV on the La 
/ states. Omitting the La / U changes the results, for 
example reducing the lattice distortions by « 50 %. For 
Sr and Ti, we use PAW potentials in which semi-core s 
states are treated as valence states, (Sr sv and Ti sv in the 
VASP distribution) while for La and O, we use standard 
potentials (La and O in the VASP distribution), and we 
use a 4 x 4 x 1 fc-point grid and an energy cutoff 500 eV. 
The lattice constants a and b are fixed to the experimen- 
tal value for cubic SrTiC>3 (3.91 A) which is the substrate 
used in the experiments |4|. The c axis lattice constant 
and atomic z coordinates are adjusted, while retaining 
the tetragonal symmetry of the crystal, until the forces 
on the ions are less than 0.01 eV/A. 

Figure^shows our calculated relaxed lattice structures 
for two representative cases: (1-8) and (2-7) heterostruc- 
tures. The largest structural relaxations occur in the 
Ti0 2 layer at the LaTi0 3 -SrTi0 3 interface (Ti . 5 in the 
upper panel and Tii in the lower panel) with the Ti dis- 
placed from its ideal position by 0.15 A in the (1-8) case 
and 0.18 A in the (2-7) structure. As a consequence, 
the lengths of the Ti-Ti distances across the LaO planes 
(Ti_o.5-Tio.5 and Tio-Ti±i) are approximately 2 % larger 
than those across SrO planes. This "ferroelectric-like" 
distortion produces a local ionic dipole moment which 
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FIG. 1: Calculated optimized lattice structures of 
superlattices [LaTi0 3 ]i [SrTi0 3 ]s (upper figure) and 
[LaTi0 3 ]2[SrTi0 3 ]7 (lower figure); half of the unit cell 
is shown in each case. The optimized z axis lattice constants 
are c = 35.09 A and 35.17 A for [LaTi0 3 ]i [SrTi0 3 ] 8 and 
[LaTi0 3 ]2 [SrTi0 3 ]7, respectively. The intertitanium dis- 
tances (lower lines) and displacements of the Ti ions relative 
to the O2 planes (upper lines) are also indicated. The center 
of the LaTi0 3 region is taken as the zero of the z coordinate 
in each case and the Ti ions are labeled by their relative z 
positions. 
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screens the Coulomb field created by the substitution of 
Sr 2+ by La 3+ ions. Moving further away from the in- 
terface, the magnitude of the ferroelectric-like distortion 
decays rapidly, while the Ti-Ti distance reverts to a con- 
stant value very close to that in bulk SrTiC>3. 

An important quantity for physical insight and more 
detailed theoretical analysis is the spatially resolved con- 
duction band charge density: loosely speaking, the Ti 
d occupancy. To obtain this, we make use of the fact 
that the ground state within LDA+J7 is a highly po- 
larized ferromagnetic state in which the magnetization 
density can be ascribed to the conduction bands. Fol- 
lowing Ref. [2^, we compute a smoothed magnetization 
density m(z), shown as the light gray line in Fig. [SJ by 
planar averaging and smoothing in z over a range ±o/2. 
We identify the integral of rn(z) over a unit cell with the 
conduction-band charge density in that cell. The total 
(summed over all cells) conduction-band charge obtained 
in this way is within ~ 1 % of the expected 1 electron per 
La ion in the two-layer heterostructure, but only ~ 0.85 
electrons per La in the one-layer structure probably be- 
cause the one-layer structure is not fully spin-polarized. 
Therefore in the latter case we renormalize the density 
appropriately. 

Fig. compares our calculated Ti d charge densi- 
ties for relaxed (black squares) and unrelaxed (white 
squares) lattices for [LaTiOaJi [SrTiC^g (upper panel) 



FIG. 2: Charge and magnetization densities for (1-8) and 
(2-7) heterostructures. Filled and open squares: conduc- 
tion band charge densities per unit cell, relaxed and un- 
relaxed superlattices, respectively, obtained as described in 
test. Light lines: smoothed magnetization densities of re- 
laxed heterostructures. Dotted, broken, and dash-dot lines: 
results of model Hartree-Fock calculation with realistic band 
parameters for unrelaxed and relaxed structures. [7 = 5 and 
J = 0.64 eV, and dielectric constant e indicated. 



and [LaTi03]2[SrTiC>3]7 (lower panel). In both superlat- 
tices, the screening provided by lattice relaxation reduces 
the charge density on the central Ti layer and produces a 
long "tail" in the charge distribution, extending far away 
from the interface. The effect is particularly large in the 
two-layer structure, reducing the middle-layer density by 
almost a factor of two. We note that in each case the 
interface layer (Tio.5 or Tii) remains electronically well 
defined, with the density dropping by approximately 0.3 
electrons between the Ti at the interface, and its neigh- 
bor surrounded by two SrO layers. The relaxed-latticc 
charge densities agree within experimental uncertainties 
with the Ti 3+ values measured by Ohtomo et al. 4] , both 
in terms of peak values (experiment: 0.3 in one-layer 
and 0.4 in two-layer) and of the slow decay away from 
the central region. 

In addition to changing the charge-density profile, as 
shown in Fig. [21 the changes in interatomic distances as- 
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sociated with the lattice relaxation lead to changes in the 
orbital overlaps. To address the effect of these changes 
in electronic structure on the many-body physics, we 
use the LDA+U results to derive a tight-binding model, 
which we solve in the Hartree-Fock approximation, com- 
puting relaxed and unrelaxed heterostructures. The form 
of the tight-binding model is discussed in Ref. ^ij. It 
has three classes of parameters: the electronic struc- 
ture, parameterized by level splittings and hoppings, 
which we obtain by fitting to our LDA+U calculations, 
"Kanamori" multiplet interactions, which we treat as ad- 
justable parameters, and the screening of the long-ranged 
Coulomb interaction, which we parameterize by a dielec- 
tric constant e chosen to approximately reproduce the 
charge densities of Fig. 

Following previous work |l4| . we represent the elec- 
tronic structure by a tight-binding model involving three 
orbitals labeled a = xy,xz,yz representing three ti g - 
symmctry Ti-0 antibonding bands, with on-site energies 
e a and nearest-neighbor hopping t s a along 5 — x,y,z di- 
rections given by the usual Slater-Koster rules, (further- 
neighbor hoppings are factors of 5-10 smaller) The 
electronic-structure Hamiltonian is H e i — ^2i[H\\(zi) + 
H±(zi)] with 

H \\( z l) = E [- 2t a( z i) cosk x ~ 2t v a (zi) cos k y + e a (zi)] 



(i) 

(2) 



a, fcii 



The heterostructures shown in Fig. ^ have too many 
bands for a direct tight-binding analysis to be practi- 
cal. We observe that there are three kinds of Ti site 
in the heterostructures (those surrounded by La, sur- 
rounded by Sr, or with two La and two Sr neighbors), 
and correspondingly three types of Ti-Ti bonds. Since 
the hopping parameters depend on the Ti environment, 
we obtain them by fitting bands obtained from LDA+U 
calculations for simplified two- layer heterostructures with 
atomic positions fixed to those found in the full (1-8 and 
2-7) structures. Hopping parameters between different 
kinds of Ti along z direction in the heterostructures are 
estimated considering the fact that the hoppings are in 
fact "second order" processes via O p state. During these 
additional band structure calculations, we eliminate the 
on-site Coulomb interaction for Ti d states which may 
cause an additional level splitting among d states, while 
the on-site Coulomb interaction for La / states is pre- 
served. 

Representative values for t s a and e a are given in ta- 
ble D We see that the largest effect of the relaxations 
is a decrease in the hoppings for d xz and d yz orbitals 
at the interface layers, and that the larger distortions in 
the two-layer heterostructure lead to larger effects. De- 



TABLE I: Tight-binding parameters for symmetry- 
inequivalent hoppings in units of eV derived by fits to 
LDA+U band calculations of simplified heterostructures 
with atomic positions taken from LDA+U calculations of 
relaxed and unrelaxed [LaTiOsJi [SrTiOsJg (labeled 1) and 
[LaTiC>3]2 [SrTiOa]7 (labeled 2) heterostructures. Parameters 
of relaxed (R) and unrelaxed (U) heterostructures are 
shown in upper and lower lines, respectively. For larger zi, 
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crease in the in-plane hoppings for d xz and d yz is as large 
as w 30 % in the (2-7) heterostructure and that in the 
z-direction hoppings is w 20 % in both the (1-8) and (2- 
7) heterostructures. Similarly, a moderate level splitting 
(~ 0.13 eV) occurs in the transition layer. 

To extract appropriate values for the dielectric con- 
stant, e, we add the long-ranged Coulomb terms to the 
tight-binding model (for details see Ref. 14]), and ad- 
just e to fit the LDA+U charge densities, as shown in 
Fig. [3 One sees that the results for the unrelaxed case 
are well described by a £ of 4, and that the effects of 
screening from the lattice relaxations can be simulated 
reasonably well (but not perfectly) by increasing the di- 
electric constant to w 15. The main deficiency is that the 
dielectric constant model overestimates the rate at which 
charge density decays far from the heterostructure, in the 
n < 0.05 region, for the (1-8) heterostructure (cf. upper 
panel). We associate this with the nonlinear screening 
characteristic of nearly ferroelectric materials. 

Finally, to investigate the changes due to the lattice 
relaxations, we compute the many-body phase diagrams 
for the relaxed and unrelaxed cases using the band pa- 
rameters derived above and the methods of Ref. 0] . To 
avoid complexity coming from the interference between 
different La regions we use isolated heterostructures in 
which n LaTiC>3 layers were sandwiched by a semi-infinite 
number of SrTi03 layers. Figure [3] compares our phase 
diagrams for relaxed and unrelaxed heterostructures. 

First we comment on features that are common to both 
relaxed and unrelaxed cases. The present calculations are 
an improvement over previously published 14] Hartree- 
Fock calculations in that they allow for the possibility 
of a fully alternating (both within plane and from plane 
to plane) antiferro-orbital ordering state. We denote this 
phase as OO-G, and see that it is in fact the ground state 
over wide regions of the phase diagram in both cases. 
This phase is favored by strong correlations (large U) and 
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FIG. 3: Spin and orbital phase diagrams of realistic three- 
band model heterostructures as a function of the intraorbital 
Coulomb interaction U with interorbital Coulomb interaction 
U' = U—2J and exchange interaction J = 0.6 eV and La layer 
thickness n computed within the Hartree-Fock approxima- 
tion. Upper panel: relaxed structure with e = 15 and appro- 
priate band parameters obtained from LDA+t/, lower panel: 
unrelaxed structure with s = 4. PM: paramagnetic, FM: 
ferromagnetic states, AF: antiferromagnetic state in which 
the magnetic moment alternates from plane to plane, OD: 
orbitally disordered state, OO-G: orbitally ordered state in 
which xz and yz orbitals alternate in x,y and z directions, 
Oxy(xz): orbitally ordered state in which xy (xz or yz) oc- 
cupancy is predominant, and 00: orbitally ordered state in 
which xz and yz orbitals alternate in z direction. 



charge densities, but with isotropic hoppings and no level 
splitting. Note that when the previous figure was con- 
structed, only calculations with in-plane translation in- 
variance were possible, the OO-G phase could not be 
studied and a higher energy layered orbital state (xz- 
yz) was found at all thicknesses. We see that the effect 
of the symmetry breaking, larger hopping intensities for 
d xy band, is to replace this phase by the ferro-orbital Oxy 
for small n. We note that the Oxy phase is not observed 
in the unrelaxed case at n = 1, even with the symmetry- 
broken hopping parameters because, in the unrelaxed 
n = 1 case, the central Ti layer (Ti±o.s) density is too 
high. Instead, the system prefers other orbitally ordered 
states with narrower xy-plane band widths. 

To summarize, we have performed first-principles cal- 
culation on LaTiOs/SrTiOs superlattices and investi- 
gated the effect of lattice relaxation on the charge screen- 
ing and change in the model parameters. Our primary 
observation is a large polar distortion of TiOg octahe- 
dra at the near La region in which Ti and O ions are 
displaced, leading to screening which reduces the central 
layer charge density, substantially consistent with exper- 
iment |4j and disfavoring staggered orbital orderings. In 
addition, for thin heterostructures at strong correlation, 
the symmetry breaking due to the distortion favors a 
novel uniform orbital ordering Oxy not found in bulk. 
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tral Facilities supported by the MRSEC Program of the 
National Science Foundation under award No. DMR05- 
20415. 



electron densities near one (24|. In thicker heterostruc- 
tures, at the edge of the high density region, we expect 
this phase to be replaced by a ferro-orbital phase for rea- 
sons similar to those discussed in Ref. |28j; this will be 
discussed in detail elsewhere [29). 

The most striking differences between the phase dia- 
grams result from the different density profiles. In par- 
ticular the lower central-layer charge density and wider 
transition regions of the relaxed structure shift the phase 
boundaries to higher U values, and disfavor the OO-G 
phase. The lattice relaxation also changes the local envi- 
ronment at each Ti site, but in the calculations presented 
here these effects cannot be isolated from the larger ef- 
fects due to changes in screening. Insight into the im- 
portance of the symmetry breaking at fixed screening is 
gained by comparing Fig. |3| to our previously published 
Fig. 1 of Ref. 14]-1 and Fig. 2 of Ref. [3-1. These 
figures show a phase diagram corresponding to similar 
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